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چکیده کامپوزیت‌های پلیمری فلزی-یونی (/(۳) مواد هوشمندی هستند که سانحتارشان شامل دو بعش فلزی و الاستومری است. از نظ رکاربردی در دسته 
پلیمرهای الکترواکتیو قرار می‌گیرند که گزینه مناسبی برای استفاده به عنوان سنسور و عملگرهای نرم با ولتاژ تحریک اندک هستند. د رکاربرد به عنوان عملکر 
خمش 1۳7/0 علاوه بر ولتاژ تحریک, متثر از عوامل ذاتی و محیطی است. پیش‌بینی #ثر پسآرامش) چالشی جدی در به کارگیری عملی 1۳0 است. اصطلاح 
اثر پ سآرامش برای کاهش تدریجی خمش 1۳1/0 و بازگشت آن به طرف کاند, تحت تحریک با ولتاژ ثابت به کار می‌رود. در این پژوهش برای اولین بار» اثشر 
پ سآرامشس در رفتار خمشی ۳(/60[ در فضای سه‌بعدی با نرمافزار 00۷501 مدل‌سازی شده است. مقایسه نتایج با داده‌ها ی آزمایش عملی, حاکی از دفت ٩۳‏ 
درصدی مدل‌سازی در پیش‌بینی جابه‌جایی است که سایر اطلاعات مهم عملکردی مانند تغییرات غلظت و توزیع تنش در ماده را نیز تأیید می‌کند. استفاده از این 
مدل برای پیش‌بینی رفتار 1۳(6 در مفایسه با انجام آزمایش‌های عملی؛ از لحاظ وقت و هزینه بسیار مفرون به صرفه خواهد بود. 
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مقد مه 

کامپوزیت فلزی پلیمری-یونی ( ۲ع۱ عصرراه۳ عنوم1 
علنو00۳0) 1۳0 از یک پلیمر یونی با گروه‌های آنیونی 
تشکیل شده است که با کاتیون‌های متحرک محلول در آب (یا 
حلال دیگر) در تعادل هستند. کامپوزیت‌های پلیمری فلزی- 
یونی (۳(۷0آها) گروهی از پلیمرهای الکترواکتیو هستند که به 
دلیل ساختار ویژه خود. گزینه مناسبی برای استفاده به‌عنوان 
سنسور و عملگرهای نرم با ولتاژ تحریک اندک هستند. یک 
نوار ]1۳1۷ از یک لایه پلیمری پوشیده شده با دو لایه رسانای 
فلزی تشکیل می‌شود. 1۳216 ماده‌ای است که حرکت مکانیکی 
را به سیگنال‌های الکتریکی (در حالت حسگری) و سیگنال 
الکتریکی را به تغییر شکل فیزیکی (در حالت عملگری) تبدیل 
می‌کند. رفتار ماده در غالب دو پدیده انتقال الکترومکانیکی و 
مکانوالکتریکی توضیح داده می‌شود. در انتقال الکترومکانیکی 
انرژی الکتریکی به مکانیکی و در انتقال مکانوالکتریکی انرژی 
مکانیکی به انرژی الکتریکی تبدیل می‌شود [1,2. 

روش پیشنهادی این مطالعه بر اساس حرکات رفت و 
برگشتی سریع مولکول‌های آب آزاد در غشای نفیون بین آند و 
کاتد است. برای درک رفتار خمشی مواد 1۳210 ابتدا یک نوار 
در فضای دوبعدی مدل‌سازی شده است تا با دریافت نکات 
کلی و جزئی از حل معادلات فیزیکی مربوط به هر ناحیه بتوان 
مدل را به فضای سه‌بعدی تعمیم داد. هدف از شبیه‌سازی نوار 
0 در اصل مشاهده تغییرات غلظت به صورت تصویر 
متحرک در نواحی مختلف ماده و نیز تجمع بار الکتریکی است 
در سطح الکترودها ایجاد خواهد شد؛ زیرا قطعا در فضای 
دوبعدی بعد سوم که مربوط به سطوح الکترودهاست قابل 
مشاهده نخواهد بود. 


ساختار پلیمر 
پلیمرهای الکترواکتیو 1۳۸۸6 از یک غشای تبادل یونی 
پرفلورینه تشکیل شده که از نظر شیمیایی با فلزی نجیب مانند 
طلا. پالادیوم» پلاتین و نقره ترکیب شده‌اند. در شکل ۱) یک 
تکار ای سول وهای بو وی مره تاد 
در ساخت این مواد هوشمند نشان داده شده است. در این 
شکل» 2 مقداری بین ۵ و۱۱ دارد و 9 تقریبا برابر ۱ (و 
شمارنده یون جنآ بهآل + و .) است. یکی از خواص 
جالب این ماده توانایی آن در جذب مقادیر زیادی حلال‌های 
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قطبی (مانند آب) است. یک پلیمر یونی پرفلورینه. ۲200 
معروف است که در اواخر دهه ۱۹۶۰ کشف شد و در اوایل 
دهه ۱۹۷۰ با فرمول شیمیایی 07۳1۳13059.03۳7) به ثبت 
رسید [1]. نفیون اولین پلیمر از دسته پلیمرهای مصنوعی با 
خواص یونی است. که به آنها آینومر می‌گویند. خواص یونی 
منحصر به فرد نفیون نتیجه ترکیب گروه‌های پرفلورووینیل‌اتر 
است که به گروه‌های سولفونات بر روی اسکلتی از 
تترافلوئورواتیلن (۲۳۲ (عممارطاهمم ۵021 /(۳۵۱) _ خانمه 
می‌یابند. نفیون به دلیل پایداری حرارتی و مکانیکی عالی» در 
غشاهای تبادل پروتون توجه زیادی را به خود جلب کرده است 
[1-3]. 


۹ ستی سس لو رتم 
"... تج تست 
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شکل ۱ ساختار مولکولی شیمیایی 2/02 (برگرفته از [1]) 


جریان یونی در پلیمر 

در حالت عملگری؛ کوج کاتیون‌های هیدراته شده در ساختار 
پلیمر که ناشی از اعمال سیگنال الکتریکی است باعث حرکت 
مکانیکی 1۳*16 می‌شود [14]. به عبارت بهتر 1۳۸66 شامل 
یک شاخه پلیمری است که به گروه‌های آنیونی و یک حلال 
دارای کانیون‌های متحرک متصل‌اند. مدل ادراکی در شکل (۲) 
نمایش داده شده است. زمانی که ولتاژی به الکترودهای 1۳160 
اعمال می‌شود بر اثر میدان الکتریکی ایجاد شده. شارش 
کاتیون‌ها (يا جریان یونی) رخ خواهد داد. در مورد 1۳۸16های 
با حلال آب کوج کاتیون‌ها ملکول‌های آب را نیز (به علت 
تغییر فشار اسمزی و تورم پلیمر در نزدیکی کاتد و فشردگی آن 
در نزدیکی آند) همراه با خود جابه‌جا می‌کند؛ که به نوبه خود 
منجر به خمیدگی ماده به سمت آند می‌شود [12]. 

در این پژوهش. شبیه‌سازی یک نوار عملگر 1۳10 
(۸۵۱۵۱0۲) ارائه شده است. به منظور تعریف مشخصات نوار 
6 در حالت عملگری شبیه‌سازی نیروهای تغییر شکل- 
دهنده و خارج صفحه‌ای بر روی نوار 1۳116 انجام شده است. 
تا کنون مجموعه‌ای از معادلات اساسی (شامل معادلات شار 
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هر دو پدیده ارائه و آنالیز شده است [1,2]. میدان الکتریکی 
اعمالی یون‌های متحرک را درون غشای پلیمری جابه‌جا می‌کند. 
در نتیجه. جریان الکترواسمزی در داخل پلیمر ظاهر می‌شود و 
متعاقبا حرکت و تغییر شکل در عملگر اتفاق می‌افتد [2]. 


ملکول آب ب 
قسمت های بدون تحرک 
(زنجیره های استخوانی) 


شکل ۲ مدل مفهومی 1۳۷6 «برگفته از [13]) 


علی‌رغم مزایا و ویژگی‌های امیدوارکننده کامپوزیت پلیمر- 
فلز پونی (1۳10) این ماده هوشمند دارای معایبی است که 
باید تا حدامکان برای بهبود عملکرد آن برطرف گردند. همان- 
طور که توضیح داده شد یکی از محدودیت‌های 1۳16 پدیده 
پس‌آرامش است که در زمان بازگشت این ماده به حالت 
استراحت اتفاق می‌افتد. پدیده بازگشت به حالت استراحت 
(حمتنیه1ع۳ 201) مشاهده شده در نوارهای با آرایش پایه 
60 که با ولتاژ 96 تحریک شده‌اند. به دلیل وجود کاتیون- 
های غیرهیدارته در شبکه پلیمری است. بنابراین» از آنجا که 
کامین‌های هش رنه از سس الکترود ایک به سمت الگیرود کاند 
در حرکت‌اند. به طور هم زمان برخی مولکول‌های آزاد آب را؛ 
به عنوان جرم اضافی به سمت کاتد منتقل می‌کنند. و اگر ولتاژ 
تحریک 90 باشد. به علت تعادل استاتیک انتظار می‌رود 
60 در نقطه انتهایی فعالیت خودش بماند و به عقب (یعنی 
حال قبل از تحریک و خمش) برنگردد. اما در واقع نوار 1۳۸66 
به عقب خمیده می‌شود و در نتیجه پدیده 1612721108 13201 را 
از خود نشان می‌دهد. پیش از این برای از بین بردن این رفتار 
راهکارهای مختلفی ارائه شده است. از جمله چندین تکنیک 
مبتنی بر استفاده از یک حلقه کنترل بازخورد برای کاهش اثر 
پس‌آرامش پیشنهاد شده است. با این حال» این تکنیک‌ها تنها 
زمانی موثر هستند که جابه‌جایی خمشی 1۳16 کوچک باشد 
[1-6. تا کنون در زمینه مدل‌سازی فیزیک پدیده انتقال 
الکترومکانیکی 1۳160 تلاش‌های قابل‌توجهی انجام شده است. 
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یک روش طبقه‌بندی مدل‌های موجود. آشنایی با نحوه توصیف 
فیک مسقله است [1], دسته اول: مدل‌هاق ارتباظدهتنده پیت 
جابه‌جایی و جریان که تقریبا به طور تجربی شکل گرفته‌اند. 
اغلب بر اساس توضیح مدار معادل الکتریکی است. دسته دوم 
مدل‌هاه شار یونی درون مواد را در نظر می‌گیرند [1-5], 
به‌طورکلی مدل‌های حالت عملگری 1۳ که تحریک 
الکتریکی خارجی را به تغییر شکل مکانیکی مواد مرتبط 
می‌سازند. می‌توان در سه گروه کلی مدل‌های جعبه‌سیاه جعبه- 
خاکستری و جعبه‌سفید دسته‌بندی نمود [1]. مدل‌های جعبه‌سیاه 
رابطه مستقیمی بین ورودی-خروجی سیستم نمایش می‌دهند که 
می‌تواند با استفاده از معادلات جبری,» اطلاعات و داده‌های 
آزمایش و بدون دانش قبلی راجع به فیزیک مسئله ایجاد شود 
[1]. مدل‌های جعبه‌سفید برای هر قسمت از مدل. مفهومی 
فراهم می‌آورد و رفتار میکروسکوپی یک سیستم را با استفاده از 
مدل‌های منحصر به زیرشاخه‌های فرایند اصلی شرح می‌دهد. به 
بیان دیگر مسئله را به مسائل ریز می‌شکند. مدل‌های جعبه‌سفید 
برای عملگر 1۳0 انتقال الکترومکانیکی را از طریق 
فرایندهای مکانیکی. الکتریکی. شیمیایی شرح می‌دهد. در 
مقایسه با سایر مدل‌هاء دستیابی به اين مدل دشوارتر است و 
برای کارپردهای مهندسی به‌ندرت استفاده می‌شود [1]. مدل‌های 
جعبه‌خاکستری بر درک فیزیک‌هایی حاکم بر رفتار سیستم 
استوارند و بخش‌های نامشخص را با استفاده از مدل‌های تجربی 
جبران می‌سازد. بنابراین این مد‌ها ترکیبی از مدل‌های 
جعبه‌سیاه و جعبه‌سفید هستند. مدل‌های جعبه‌خاکستری 
به‌سادگی از نقطه‌نظر شناسایی سیستم یک مصالحه بین نیازهای 
اساسی مدل برای درک سیستم برقرار می‌کند [1]. با توجه به 
اینکه شبیه‌سازی رفتار 1۳۷ نیازمند مطالعه دقیق و یافتن 
ارتباطی بین خواص مکانیکی و شیمیایی فلز و پلیمر در صورت 
اعمال میدان الکتریکی است. تمرکز اين مقاله روی مدل‌های 
تحلیلی (از انواع مدل‌های جعبه‌سفید) [1] است. هم‌چنین چون 
تفاوت جنس خواص مکانوالکتریکی و شیمیایی در الکترودها 
و نفیون و شرایط مرزی بین دو ماده نیازمند حل میکروسکوپی 
معادلات اساسی مشترک است. از میان انواع مدل‌های تحلیلی به 
مدل‌های فیزیکی و شبیه‌سازی رفتار 1۳۸60 با استفاده از 
نرم‌افزار کامسول که اساس حلگرهای (50176) آن بر پایه 
روش‌های المان محدود (۱6/0006 تصصعاظ ننهز۳) است. 
می‌پردازد [1.3]. نرم‌افزار چندفیزیکی کامسول یک مجموعه 
کامل شییه‌سازی. است که می‌تواند مععادلات دیفرانسیل 


سیستم‌های غیرخطی را توسط مشتق‌های جزئی به روش اجزای 
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محدود (۳/۷) در فضاهای یک دو و سه‌بعدی حل نماید. اين 
نرم‌افزار می‌تواند در حضور چالش‌هایی نظیر میدان‌های 
الکترومغناطیسی. کشش. دینامیک سیالات و دینامیک گاز به 
خوبی راهگشا باشد. کامسول همچنین فرصت حل مسئله با 
فرمول‌های ریاضی (در فرم معادلات) و فیزیک (انتخاب مدل 
فیزیکی) را به کاربر می‌دهد. بدیهی است. در هر مسئله یک 
سیستم معادلات حل خواهد شد اما تفاوت فقط در نحوه 
استفاده از سیستم‌های فیزیکی نهفته‌است. می‌توان از معادلات از 
پیش تعریف شده برای اکثر پدیده‌ها در علوم و فناوری, مانند 
انواع روش‌های انتقال حرارت. برق» تثوری الاستیسیته. نفوذ و 
انتشار مولکولی. انتقال جرم. انتشار موج و جریان سیال استفاده 
کرد. 

در ادامه بنا بر آنچه گفته شد. ابتدا در بخش روش تحقیق 
مروری بر معادلات حاکم بر فیزیک 1۳16 و نحوه مدل‌سازی 
کامپوزیت در حالت سه‌بعدی را با جزئیات بررسی خواهیم کرد. 
سپس نتایج حاصل از شبیه‌سازی را به صورت گرافیکی در 
قالب نرخ تغییرات غلظت. تنش, نیروهای بدنه. تغییرات 
جابه‌جایی نسبت به زمان» تحت تحریک ولتاژی بررسی کرده و 
از نظر کمی نتایج حاصل از مقایسه جابه‌جایی نوک 1۳10 را با 


داده‌های ثبت شده از آزمایش عملی ارزیابی می‌نماييم. 


اثر الکتریکی الکترودها 
مطالعات انجام شده حاکی از آن است که جریان الکترودها 
روی رفتار انتقال الکترومکانیکی تأثیررگذار است و الکترودها به 
طور قابل توجهی بر دینامیک بار و عملکرد عملگری 1۳316 
تأثیر می‌گذارند. برای در نظر گرفتن الکترود در مدل. جریان 
یونی در پلیمر با جریان الکتریکی در الکترودها کوپل می‌شود. 
برخلاف مدل‌های فیزیکی پیشین. الکترودها ایدئال نیستند اما از 
هدایت الکتریکی با ضریب انتقال محدود برخوردار هستند. 
اگرچه اثر الکترود در مورد انتقال الکترومکانیکی (مقاوم در برابر 
جریان) بیشتر قابل درک است. اما فیزیک پایه برای هر دو نوع 
انتقال یکسان است. معادله (۱) جریان یونی را در بخش 
پلیمری 1۳260 توصیف می‌کند. در مورد الکترودهاه شکل 
دیفرانسیلی قانون اهم برای چگالی جریان به صورت رابطه (۱) 
است [2]: 
)۱ [- < 6۸۷ 
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که در آن ۵ هدایت الکتريٍ کی الکتر ها ا ست. ۷ و [ به 
ترتیب پتانسیل الکتریکی و چگالی جر یان در الکترود ها | ست 
[1 
عوامل موْثر بر رفتار 1۳۸16 
به طور کلی عوامل مختلفی رفتار 1۳۷6 را در هر دو حالت 
سنسوری و عملگری تحت تأثیر قرار می‌دهند که می‌توان آنها را 
در دو دسته کلی عوامل ذاتی و محیطی قرار داد. تأثیرگذاری 
همین عوامل است که سبب می‌شود مدل‌های خطی انعطاف- 
پذیری کافی برای مدل‌سازی رفتار 1۳*16 را نداشته باشند و 
رفتار ورودی و خروجی در این پلیمر تنها با استفاده از مدل‌های 
غیرخطی قابل مدل‌سازی و پیش‌بینی باشد. آزمایش‌های انجام 
شده بر روی 1۳1۷0 نشان می‌دهد. میزان خمش در حالت 
عملگری و درک جابه‌جایی سنسوری این ماده تحت تأثیر 
عواملی قرار دارد که متناسب با کاربرد مورد نظر در ساخت این 
ماده باید در نظر گرفته شود. از اين پارامترهای موثر که عمدتا 
متأثر از نحوه ساخحت 1۳10 هستند. اصطلاحا تحت عنوان 
عوامل ذاتی نام یاد می‌کنیم. مهم‌ترین عوامل ذاتی عبارتند از: 
. نوع یون‌های لایه پلیمری (غلظت بالاتر من جر به پا یداری 
خمش می‌شود). 
۲ ضحخامت یون‌های لایه پلیمری (هرجه ضخامت بیشتر 
انعطاف و خمش بیشتر). 
۳ ابعاد نوار 1۳116 (نسبت طول به عرض. طول بٍ شتر سبب 
خمش بیشتر می‌شود ولی برای عرض برعکس است). 
۴ جنس الکترودهای پوشاننده لایه پلیمری و ضخامت آنها. 
۵ فاکتور خحستگی (ویژگی‌ای از ماده که باعث می‌شود 
جابه‌جایی به مرور کمتر شود). 
تما نله 


ویژگی‌ها و کاربردهای ۲۳۸۷۲ 

6 ها به دلیل سبکی, انعطاف‌پذیری نرمی, سازگاری با 
محیط یونی داخل بدن و شباهت رفتاری با رفتار عضلات بدن 
از جمله مواد پرکاربرد و مورد مطالعه در علوم مهندسی و 
خصوصا مهندسی پزشکی به شمار می‌روند و تحت عنوان 
باهتخههای عیاقو 111 

به علت انعطاف‌پذیری بسیار خوب و اعمال نیروی نسبتا 
زیاد (نیرویی نزدیک به ۴۰ برابر وزن خود) این مواد در مقابل 
اعمال ولتاژ بسیار کم و عکس‌العمل بسیار سریع آنها (در کسری 


یشیرق راز زار هه ۱۴۶۲ 


آرزو حراجی, نادیا نقوی, حجت زامیاد 


از ثانیه) در ساخت ماهیچه‌های مصنوعی کاربرد دارند. مطابق 
توضیحات ارائه شده ویژگی‌های منحصر به فرد 1۳1۷6 باعث 
شده در کاربردهای گوناگون صنعتی و پزشکی مورد توجه باشد 
که به برخی از آنها اشاره شد. در این قسمت به طور خحاص 
مهم‌ترین کاربردهای پزشکی به تفکیک مورد بررسی قرار می- 
گیرد [27,28]. 


کاربردهای پزشکی 1۳16 در حالت عملگر 

6 به دلیل انطباق زیاد. ولتاژ فعال‌سازی پایین 
توانایی کار در شرایط محیطی مختلف و سهولت ساخت و 
سفارشی سازی. نویدبخش فعال‌سازی. سنجش و برداشت 
انرژی در رباتیک نرم. مهندسی زیست پزشکی, ذخیره انرژی و 
حتی کاربردهای فضایی هستند [3-27]. 

از مهم‌ترین دغدغه‌های استفاده از ابزارآلات پزشکی 
محرک در تشخیص و درمان بیماری‌های داخلی بدن انسان 
سازگاری با محیط الکترولیتی بدن, تأمين انرژی لازم و ابعاد 
قابل تغییر برای ایجاد حرکت دل‌خواه در عنصر محرک مورد 
نظر است. همین مسئله باعث شده 0آها به دلیل تطابق با 
محیط الکترولیتی بدن, ولتاژ کاری کم و امکان استفاده با ابعاد 
مختلف یکی از بهترین گزینه‌ها به عنوان یک اکچویتر باشند. 
۱. ماهیچه‌های کمک‌رسان قلبی (بطن مصنوعی) 
۲. عضلات صاف و اسفنکتری مصنوعی 
۳. اصلاح عیوب انکساری چشم 
۴ پروتزهای مایو الکتریک دست 
۵ استفاده به عنوان دریچه دستگاه تحویل هدفمند دارو به 

صورت بی‌ سیم 

۶ انگشت مصنوعی [27]. 


بیان مسئله اثر خستگی در 1۳[۷۲6ها 
پدیده بازگشت به حالت استراحت (1612:۵1:08 0۲ظ) 
مشاهده شده در تير 1۳*16 که با اعمال ولتاژ 96 تحریک 
شده‌اند. به دلیل وجود ملکول‌های غیرهیدارته در شبکه پلیمری 
است. بنابراین از آنجا که کاتیون‌های هیدراته از سمت الکترود 
آند به سمت الکترود کاتد در حرکت‌اند. به طور هم‌زمان برحی 
مولکول‌های آزاد آب راء به عنوان جرم اضافی به سمت کاتد 
منتقل می‌کنند. و اگر ولتاژ تحریک 0 باشد. به علت تعادل 
استاتیک انتظار می‌رود 1۳10 در نقطه انتهایی فعالیت خودش 


ال سای یار فجاره صمف ۱۳۰۲ 
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بماند و به عقب (یعنی حال قبل از تحریک و حمش) برنگردد. 
اما در واقع نوار 1۳146 به عقب خمیده می‌شود و در نتیجه 
پدیده 161221108 3201 را از خود نشان می‌دهد. حال تصمیم 
بر آن است که ابتدا برای درک رفتار حمشی این مواد یک نوار 
6 را مشابه [2] مدل‌سازی نموده و با دریافت نکات جزئی 
سل ابعادلات فیزیکن موبوط به هر تایه تواشیم مدل ,زا به 
فضای سه‌بعدی تعمیم دهیم. 


روش تحقیق 

معادلات فیزیکی -شیمیایی 
بر حلاف انتقال مکانوالکتریک در انتقال الکترومکانیک. تحت 
تأثیر اعمال اختلاف ولتاژ به الکترودها 1۳10 دچار تغیبر شکل 
و بخش آزاد تیر از انتها شروع به خمیدگی می‌کند. آنچه حائز 
اهمیت است این است که علت اصلی هر دو پدیده 
مکانوالکتریکی و الکترومکانیکی جریان یونی القاشده و در 
نتیجه توزیع بار غیرمتقارن در مجاورت الکترودهاست. جریان 
یونی در پلیمر از معادله نرنست‌پلانک طبق رابطه (۲) محاسبه 
می‌شبود 01,21 


06 
2 + ۲۰/۲6 - 2۷۳6۲۵ - ۷6۵۷۵۴( - 0 (۳) 


که در آن 0 غلظت کاتیون؛ ۸ قابلیت تحرک (ونانم]0۷) 
کاتیون. ثابت انتشار» ۳ ثابت فارادی. ,7 عدد بار ۷ حجم 
مولی (که معیاری برای کمیت‌سنجی آب‌دوستی کاتیون‌ها به 
شمار می‌رود). ط فشار حلال وه پتانسیل الکتریکی در پلیمر 
است. لازم به ذکر است که پتانسیل الکتریکی 4 در داخل پلیمر 
و پتانسیل الکتریکی ۷ در الکترودها دو متغیر متفاوتند. ضریب 
پویایی ‏ را می‌توان به‌صورت رابطه (۳) بیان کرد [1,2]: 


_ 
۲ (۳ 


در این رابطه 8 ثابت گازها و ۲ دمای مطلق است. معادله 
(۲) رابطه اصلی حاکم بر توضیح پدیده انتقال در ماده 1۳16 
است. این رابطه علاوه بر زمان تحت تأثیر گرادیان‌های میدانی 
نظیر گرادیان پتانسیل الکتریکی ۰۷4 گرادیان غلظت )۷ و 
گرادیان فشار حلال ۸۵۳ نیز است. این گرادیان‌ها در هر دو نوع 
انتقال وجود دارند. گرادیان پتانسیل ۵ با معادله پوآسون (رابطه 
(۳)) توضیح داده می‌شود. در روابط (۴) و (۵) م چگالی شارژ 
غلظت آنیون و متغیر 8 ثابت گذردهی مطلق موثر است که 


قی‌توان: انوا به: ضورت. ۳68 پوفنت. که دور آ۵ و6 
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ضریب گذردهی خلاً و برابر مقدار (/۸/۸۵۱۰۱۲۳ است. 
که در آن مدا غلظت اولیه کاتیون با آنیون است [1,2]: 


)۳( ل - ۳2۵ 
: ۷6 


(() (م) )۳ < م 
غلظت کاتیون از رابطه (۶) به دست میآ ید در حالی که 
غلظت آنیون وابسته به کشش حجمی محلی است [1,2]. 


(ع ۷ 


<< 
در رابطه (۶) ۷ بردار جابه‌جایی محلی است. مقدار مثبت 
کشش حجم محلی نشانگر افزايش حجم محلی و منفی بودن 
آن به معنای کاهش حجم است. تغییرات حجم در ماتریس 
پلیمر غلظت آنیون‌ها را به عنوان بخش بدون تحرک پلیمر 
تحت تأثیر قرار می‌دهد. از اين رو متغیر غلظت آنیون ,0 با 
رابطه (۷) تعریف می‌شود و م) غلظت اولیه کاتیون/آنیون است 
[1,2]. 


0 0-1۷( 2 


برای یک 1۳*0 با ضخامت ۰/۴۵ میلی‌متر با ضریب 
گذردهی مژثر دی الکتریک ۰/۲۲۵۳/۲و غلظت اولیه کاتیون 
0/3 ۱۲۰۰ است. افت پتانسیل ۱ ولت در اين ناحیه باریک 
منجر به ایجاد گرادیان تقریبا برابر با ۱۰۶ ولت بر متر می‌شود. با 
در نظر گرفتن کشش 1۳0 خم شده در حدود ۰۱/ با مدول 
یانگ ۲۰۰ مگاپاسکال و نسبت پواسون ۰/۴۹ فشار را می‌توان 
۵ مگاپاسکال برآورد کرد. این تخمین گرادیان فشار به 
صورت قراردادی در محدوده ۱۰۶ ۲ ۳/۵ مطرح می‌شود. اگر 
عبارات ۳۲۷4 و ۸۷۷۳ در رابطه (۱) مقایسه شوند. با توجه به 
اینکه ثابت فارادی و ۸۷ به ترتیب برابر مقادیر (۹۶۴۸۵)0/01 
و (۶*۶-۱۰)۳:3/001 است. می‌توان مشاهده کرد که گرادیان 


)۸( ۱۸۷۷۲۳۱ << (۳۲۷۵| 
بنابراین» بخش تتبوم شار شکل م کانیکی ۱92۵۵ در 
محاسبات انتقال الکترومکانیکی در نظر گرفته نمی شود. علاوه 


بر اين. همان طور که توضیح داده شد. چگالی بار با یرو های 


بدنه در ماده مرتبط است. یعنی [1,2]: 
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)0( (2)0 - (0) - ۲ 
بنابراین فرآیند کوپلینگ, فعل و انف عال بین ب خش بدون 
تحرک و حلال را توصیف می‌کند و عبارت گرادیان فشار قا بل 


چشم‌پوشی است. پس معادله نرن ست‌پلا نک در مورد عم گر 
الکترومکانیکی به صورت زیر استفاده می‌شود [1,2]: 


06 1 
2 + ۲۷۰/۲6 -7۳6۲۲( < 0 ۰( 


لازم به ذکر است که برای اکثر محاسبات عملی تقریب 
و) < م) کاملا منطقی است. فقط در موارد انتقال 
مکانوالکتریکی, که ترم‌های مربوط به شار بسیار کوچک هستند 
باید از رابطه (۷) استفاده شود. این ساده‌سازی‌ها مفاهیم 
معناداری دارند. یعنی محاسبه سیستم ۳(۷۳ (پوآسون نرنست 
پلانک) از مدل الاستیک خطی جدا می‌شود. 


پیاده سازی مدل 1۳۱۷۲ در سه‌بعد 
مدل انتقال الکترومکانیکی با روش ۲۳۸۲۲ ( احعصهاط انمز۳ 
6 حل شده است. به‌منظور محاسبه فیزیک‌های مدل. 
سیستم معادلات پوآسون» نرنست‌پلانک و فرم دیفرانسیلی قانون 
اهم به‌طور هم‌زمان حل شده است. برای انجام این کار» شرایط 
مرزی مناسب تعریف شده و مدل با نرم‌افزار المان محدود 
کامسول پیاده‌سازی می‌گردد. معادلات ۳۱۴ در ناحیه پلیمری 
حل شده و چگالی جریان الکتریکی در نواحی الکترودی پیاده- 
سازی شده‌اند. حل این معادلات با انتخاب فیزیک‌های مربوطه 
انجام خواهد گرفت [1,2]. 
از آنجا که در نسخه‌های جدید نرم‌افزار کامسول امکال 

شکل‌دهی حجمی (۲۳:0۷06 مدل دوبعدی به فضای سه‌بعدی 
در اختیار کاربر قرار گرفته است. پس از بررسی حالت دوبعدی؛ 
برای شبیه‌سازی نوار 1۳10 در فضای سه‌بعدی دو روش پیش- 
روست: 
۱ از ابتدا در فضای سه‌بعدی. تعریف هندسه نوار 1۳۷ و 

سپس تحلیل معادلات مربوطه انجام شود. 
۲. پس از آنکه مدل دوبعدی به نتیجه رسید. با گسترش دادن یا 

به نوعی حجم بخشیدن به هندسه دوبعدی تعریف شده. 


مدل به فضای سهبعدی انتقال پابد و تحلیل شود. 


ایس ری هازی قیرط متفه ۱۳:۲ 


آرزو حراجی, نادیا نقوی, حجت زامیاد 


از این دو روش روش اول به دلیل امکان مشاهده سطوح 
الکترودها و آماده‌سازی هندسه برای محاسبات المان محدود به 
منظور انجام شبیه‌سازی انتخاب می‌شود. 

مدل برا تبر یک سر درگیر تهیه و بررسی شده است. 
هندسه مدل ایجاد شده در ابعاد ۲۸/۵ ۰/۴۵ و هم 
طراحی و پیاده‌سازی گردیده است. مقادیر ابعاد تير با هدف 
اعتبارسنجی نتایج شبیه‌سازی سه‌بعدی با داده‌های آزمایش 
عملی. مشابه نمونه نقوی و همکارانش در [1] در نظر گرفته 
شده اتیتتا: در این هندسه 171۳0 ۲/۵ از طول تیر حهت اعمال 
ولتاژ در یک کیره ثابت در نظر گرفته شده است, در نواحی 
الکترودی از فیزیک عاصعتت) 160۱710 برای اعمال ولتاژ 
تحریک خارجی و اثر الکتریکی الکترودها و در نواحی (بخش 
داحل گیره و بخش آزاد) پلیمری از فیزیک 08 10050071 
695 12110160 برای پیاده‌سازی معادلات فیزیکی-شیمیایی 
استفاده شده است [27]. هم‌چنین فیزیک ۱/6۵21 5010 
اعمال نیروهای 
الکترواستاتیکی و فیزیک ۳۲8 0۶ ۲0۲۳۲ ادتعصع0 برای حل 
معادلات دیفرانسیل مورد استفاده در مدل به کار گرفته شده 
است [13]. برای حل معادلات دیفرانسیل مشتقات جزئی در 
نرم‌افزار مش‌بندی پا ایجاد یک طرح شبکه‌ای در هندسه ماده 


جهت محاسبه جابه‌جایی ناشی از 


ضروری است. در دستیابی به نتایج منطقی و صحیح انتخاب 
مش و نیز عدم وابستگی نتایج به تعداد المان‌ها و نوع مش 
بندی (اصطلاحا رسیدن به مدل مستقل از مش) بسیار حائز 
اهمیت است [1]. انتخاب مش به میزان زیادی وابسته به آن 
است که کوپل (جفت شدن) جریان الکتریکی/شاریونی بر 
نتاس یه زاو شا کی سکب فان ارس آ دا برا ه 
اول به انتگرال‌گیری روی کل ناحیه نیاز است که در نرم‌افزار 
چند فیزیکی کامسول می‌توان آن را فقط در مش‌های مثلثی 
محاسبه کرد. در مورد کوپل مبتنی بر قانون گاوس مش‌های 
نگاشت (0۷1۵0108» مورد استفاده قرار می‌گیرد. مزیت مش 
نگاشت در مقایسه با مش مثلثی» کنترل بیشتر بر اندازه و توزیع 
عناصر است. صرف نظر از انتخاب. مش در نزدیکی مرزهای 
الکترود و نفیون به علت تغییر زیاد غلظت ۷ و گرادیان 
پتانسیل 7۲0 باید مش‌بندی مناسب باشد. نمونه مش‌های مثلثی 
ق تکاشتت ون شک (۳) نمایقن ذافه فده است: [1327]. مش 


های نگاشت به صورت مکعب و مش‌های دو بعدی مثلثی در 


سل مبيی ق جهن جنماره شهب ۱۴۰۲ 


۸۷ 


حالت سه‌بعدی به صورت هرم چهاروجهی (161072060721) 
فی سود: 

به‌منظور انجام حل معادلات روی فیزیک‌های مدل از دو 
مطالعه وابسته به زمان به‌گونه‌ای که بخش اول تحلیل اثر 
الکتریکی الکترودهاء فرایندهای الکتروشیمیایی و بضش دوم 
به‌منظور تحلیل معادلات ارتعاش مکانیکی. اعمال نیروی 
وابسته به زمان استفاده شده است. در این مقاله مطابق با [2] 
ثابت انتشار رل ۰۱۱-۱۰ ۰/۷ ضریب گذردهی الکتریکی 9 
۳ مدول یانگ ۶ ۴۱۷۳۵ و ضریب پوآسون ۰/۴۹ در 
نظر گرفته شتده: انستا: تمامی پارامترها و مشخصات ماده 1۳۷۲ 


به استناد [1] مقداردهی شده‌اند. 


(الف) 


کر کر ۳۹ 
و 
در ۳ 
گر 


۳ 

> 

7 بد 
هراکه ۵۲ 


کی 


لب) 


شکل ۳ نمونه‌هایی از مش‌بندی در نرم‌افزار کامسول برگرفته از [13] الف) 


مش نگاشت و ب) مش مثائی 


نشریه مهندسی متالورژی و مواد 


۸ 


نتایج 

با توجه به آنکه کوچک کردن ابعاد و افزایش تعداد المان‌های 
مش‌بندی بسته به مشخحصات پردازنده و حافظه رایانه امکان‌پذیر 
است. ملاک قابل پذیرش بودن مش‌بندی با توجه به میزان 
پارامترهای آماری نظیر 7القتاو احعصهاه ععه۸۷۵ قابل ارزیابی 
است. این پارامتر ویژگی بهینه بودن مش را با نمایش مقداری 
بین ۰ و ۱ تبیین می‌کند. این مقدار هرچه به ۱ نزدیک‌تر باشد» 
به مش ایدئال نزدیی‌تر است. مش‌بندی بهینه با استفاده از 
مجاور, تأمین حد قابل قبولی از فاکتورهای آماری تنظیم مش 
کاهش ابعاد مش لزوما نشان‌دهنده کیفیت مطلوب مش‌بندی و 
افزایش دقت محاسبات نیست. بلکه ملاک تشکیل مشی است 
که بتوان به نحو مطلوب مدل را به مقدار بهینه نزدیک و 
نزدیک‌تر کرد و از واگرایی تخمین در روش های گسسته‌سازی 
ریاضی جلوگیری نمود. 


ی‌ااعناواه 


جاجهه ع21ام۵) 
3264 جع ۷۶۱۲ ۶ع۱۷ 


عاصعاع الما عمط ۲ع۳۲ع۴۱ 
2250 ارو 
2435 9 


8 :2۳۱05 ۴۵9 
4 :۶۱2۳۱605 ۱۷۶۳۲۵۲ 
یااوناهاز ۱2۳۱۵۳۲ 0۵۳۴۵۵۱۲ 
2250 اعاع حعاع ۵۲ ۲عاطبلا 
0 ازادیان خطععصعاع بطم 
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000709 
۲ 6413 


۴۶۱۶۳۱۶۱۴ ۷۵۱۷۱۱۵۵۵ 


۱۷۶۳۱ ۵۸ 


شکل ۴ ارزیابی مش‌بندی 

تير یک سر درگیر شبیه‌سازی شده در این پژوهش. د حت 
اعمال ولتاژ مستقیم با دامنه ۵ ولت و به مدت زمان ۱۴۰ ثانبه 
تحریک شل‌ه انستت: 

در شکل (۵) به وضوح تغییرات مژلفه‌های نیروی وارد بر 
مرز میان پلیمر و الکترود مشاهده می‌گردد. در مدل دو دسته 
نیروی مرزی دیده می‌شود. پیکان‌ها نشانگر نیروهای مرزی 
ناشی از اعمال نیرو الکترواستاتیک در مرز الکترود بالایی و 


نشریه مهندسی متالورژی و مواد 


مدل‌سازی سه‌بعدی اثر خحستگی در کامپوزیت‌های پلیمری فلزی پونی 


پلیمر است و مخروط‌ها نمایانگر نیروهای مرزی ناشی از میدان 
الکترواستاتیک به حدفاصل الکترود پایینی و پلیمر است. 

با توجه به طیف رنگ گرادیانی در شکل (۶) می‌توان 
دریافت که بیشترین تنش وارد به تير 1۳1 در نواحی 
قرمزرنگ نزدیک به گیره و کمترین تنش به نواحی آبی رنگ 
(نوک تیر) اعمال می‌گردد. از ناحیه مهارشده در گیره به دلیل 
ثابت بودن صرف نظر شده است. 

شکل (۸ نتیجه تغییرات غلظت نفیون پس از تحریک 
60 را نشان می‌دهد. همان طور که در بزرگ‌نمایی این شکل 
دیده می‌شود در نمودار تغییرات غلظت نواحی الکترودی بی- 
رنگ هستند. با توجه به این تغییرات غلظت. تورم 1۳66 در 
ناحیه قرمزرنگ این شکل رخ می‌دهد. لازم به‌ذکر است 
گرادیان غلظتی دقیقا معکوس گرادیان پنانسیل در تواحی 
پلیمری است. 

مهم‌ترین نتیجه این مدل‌سازی» تحلیل جابه‌جایی نوک تیر 
60 است. در شکل (۸) داده‌های حاصل از اعمال فیزیک- 
های مختلف روی جابه‌جایی تیر 1۳۸16 مدل‌سازی شده و 
نمونه واقعی [21] با یکدیگر مورد مقایسه قرار گرفته‌اند. انطباق 
داده‌های حاصل از مدل‌سازی با نمونه مشابه واقعی تحت ولتاژ 


مستفیم ۳/۵ ولت نشان‌دهنده موفقیت‌آمیز بودن مدل طراحی 
شده ابیت بر مبنای معیار ذکر شده در رابطه (۱۱ بین داده‌های 
بیل‌ساای: 1۳0 و تمه واقعی انطیاق ٩۳‏ فرصی: وه 


دارد [1]. 


(1 - 100 - (96)ودنااز۳ 


در این رابطه 4 معرف جابه‌جایی حاصل از مدل‌سازی, 4 
جابه‌جایی نمونه واقعی. ز گام‌زمانی و 2 تعداد کل گام‌های 


زمانی ثبت شده است. 


سالن.سی ق چهار: ساره سه ۱۴۰۲ 


آرزو حراجی, نادیا تقوی, حجت زامیاد ۸٩‏ 


(0تامع) ۱۵05 800003۲۷ 4 ۶ ۶ 711۳۳6213 


(تحت اعمال ولتاژ ۲6 با دامنه ۳/۵ ولت در مدت زمان ۱۴۰ تانیه) 


(۱۱/۳۵۶) ععه۱۳: ۲/5۵5 ۷۵۵ :5۵۲6۵6۵ 4 شیر 5 11۳6-13 


شکل ۶ توزیع تنش در 1۳1۷16 . (تحت تحریک الکتریکی مستقیم با دامنه ۵ ولت در مدت زمان ۱۴۰ ثانیه) 


11۳6-13 ۶ ۳۳ ۱۳0/۳3 


0 


(1 


1 


شکل ۷ گرادیان غلظت در نواحی پلیمری کامپوزیت 1۳۷16 (تحت تحریک الکتریکی مستقیم با دامنه ۳/۵ ولت در مدت زمان ۱۴۰ ثانیه) 


سال سی و چهار» شماره سهء ۱۴۳۰۲ نشریه مهندسی متالورژی و مواد 


۹۰ 
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۱ (۸5 94 1۹ م 1 10 اصعسی‌هام‌ونط_ 
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شکل ۸ مقایسه جابه‌جایی نوک 1۳16 در حالت پس‌آرامش (تحت اعمال 
ولتاژ ۲90 با دامنه ۳۸۵ ولت در مدت زمان ۱۴۰ ثانیه) 


با نمونه مشابه مدل‌سازی‌شده در نرم‌افزار جند فیزیکی کامسول 


نتیجه گیری 
پدیده بازگشت به حالت استراحت (همناق62 01ظ) 
مشاهده شده در نوارهای با آرایش پایه 1۳۷66 که با اعمال ولتاژ 
0 تحریک شده‌اند. به دلیل وجود مولکول‌های غیرهیدارته در 
شبکه پلیمری است. بنابراین. از آنجا که کاتیون‌های هیدراته از 
سمت. الکفرود آند به. سمت الکترود: کانك در حر کنت‌انل به‌ظوز 
هم‌زمان برحی مولکول‌های آزاد آب را به‌عنوان جرم اضافیء به 


شفت 


مدل‌سازی سه‌بعدی اثر خحستگی در کامپوزیت‌های پلیمری فلزی یونی 


سمت: کاند متفل می‌کننده و اگر ولناز تحریک: 06 باشد: به 
علت تعادل استاتیک انتظار می‌رود 1۳۸0 در نقطه انتهایی 
فعالیت خودش بماند و به عقب (یعنی حالت قبل از تحریک و 
خمش) برنگردد اما در واقع نوار 1۳16 به عقب خمیده 
می‌شود و در نتیجه پدیده 161220100 13201 را از خود نشان 
می‌دهد. هدف این مدل‌سازی پیش‌بینی میزان بازگشت 1۳۷6 
به حالت پس‌آرامش بود. با توجه به نتایج گزارش شده در 
بخش چهارم و انطباق قابل قبول مدل با واقعیت. می‌توان مدل 
0 سهبعدی را به هندسه‌های دیگر با ابعاد الکترود و 
ضخامت پلیمر گوناگون نیز تعمیم داد. مدل توصیف‌شده قادر 
است علاوه بر جابه‌جایی ملفه‌های دیگری نظیر تغییرات 
غلظت. تغییرات پتانسیل» تغییرات نیروهای مرزی» و میزان تتش 
وارد بر هر نقطه از هندسه انواع 1۳146 را با در دست داشتن 
خواص مکانیکی, شیمیایی و الکتریکی بخش پلیمری پیش‌بینی 
نماید. قطعا با توجه به هزینه‌های تهیه یک تير 1۳1۷16 استفاده از 
این مدل از لحاظ وقت و هزینه بسیار مقرون به صرفه خواهد 
بود و مانع از هدررفت نمونه‌های سالم و نیز تهیه سخت‌افزار 


دقیق ثبت داده از این ماده خواهد شد. 


تقدیر و تشکر 
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تشر مچدیسی رای و مراد 
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لمامصه تمصراوم متصمد گه صمتاهتتاصم10 مسجت همه صمتمطهماع صمتصتتمولظ رامطعدل ۱۲۰ 0صه تصفاوه‌طحصصم 10[۳] 
4 ,582-598 .00 .5 .10 ,26 ۷۵۱۰ وگ ره ۷۵۹۵ ۲۱/۵۱۲ ۰۷۰ ,5تماحنتامه عالوومجیع 

حصحط م6 عوجر ۷21102610 عم تممنوآم ۳1۵2061601710 راقل1 ۲۱۱۵-۲۷ ۱۷۲ صه بصایه۷-معطمصرهه .1 بهاممع۱]20-۳:ظ .[[11] 
0۰ ,۱۷27 ,616600 .0 ,6166 ۷۵۱۰ ,۳ ۳۵6( هو اجمعتا اما ۵۲۹۱۷۲ 1-9۷۷ )0طامت0تمنجط 5120 

۲ 5008018 مالوممممه قامصه تمصصروامی عمتجم تفگ ممتاهکممممدمه 0100۷9عم 0عقعها ۵01ص ۸ ص12 ۷۲۰ 24 هدومن بع[12] 
۰ ۸0۲ ,43-49 .00 ,224 ,۷۵۱ ,.عجاظ یگ ۸0۵۵ ععی " رتصعجطمعیاعجمطط کصمجطرهمه۱ مرول 

,منمها۵و .بآ ) عک ممصحتتانظ ما بهف2700ظ نک بل۸]022 یک مفع۷ 5,008 ,0۰92002700 ۱۷۲۰ ,ما20 .0۰ :13[۸۵] 
اصععنه ه مقجمموع۲ »1۵۵0۵20 1۵۲ ممیود مصتطم‌فصه ملع ۲ ۵/۳) مه وم 0عوده وعمتهنتامه ۱۳۱۷۲ ۶0۲ ععلا0تنجمم 
2022(۰) .25 بمت2ععمک ولهتنماه۷ .01065 تصصبط ۲۵12۷۵ 

عنصم1 ممصممصم۳۵ م۲ ه ۱۵0۳۷2۵۲۵0 ۴ پطاومصمصصله۳ ,۱ کی ماو تاتممصصعط .۲ بهع7ه9 ۲۰ ۷ بتصدهاوعطامصهه .]14[1۷] 
ر2015تااعظ محصم۱ ۳۲۲۵۵-صم)۲۵ ۵۶ ممتاف‌«دامک ا0ظ ۶ه ممتامصتصصتاط عم ۱۷۲۵0۵۵0 060ته-عمصهها )ورن خر :عماقتتاع۸ 

1 ممنطممه9۵0 ,2۲ ,مصعتنایه با م210082ظ .5 ۸۲022 بک وقع7 002 با ,820602700 0۰ ,۷ رماافتاوتا7 .0 م[15] 
کصععنه 2 مقجمموع۲ ۵۵0۵20 1۵۲ ممول مصتطم‌فصه مصلعت 0/۳۲) مه وم عقدها وعمتهتتامه ۱۳۱۷۲۵ ۶0۲ ععاامتنجمم 
2022(۰) 1۵8681600,25۰ ۱۷۲2۲۵۲1215 ۰ رو01016 تحصیاط ۲۵۱۵۷۵ 

مط وه حمتاهعنامو2 معه)۷۵1 مه فصمتنومه امصمصصما هه ۵ ممصمتالما * رواواه‌طامو. ,۷۲ ۲۰ بتاو .1 بامصتم ,۷ .ه[16] 
ب0.11ظ 31 ۷۵۱۰ اه وه واهه۷/۵ ۸۵ روتمتهنتامج 1۳۸۸۵ ۲ -صمنلداط ۵۶ معصمصمممم لهمنصهههعصهمتامعام 
جشضلا2 ,112031 

00۵/۵۵۵۵۵0۵ ۵0وومن ند ه ۵۶ ممتاهم۳20 * .5 پتا۷ لگ وبا ولا معصعط2 .2 رعصقطا2 2 رعصه۷۷ ,۱۷۲ ,عصعط2 ,17[۳۲] 
,356 همعط :ظ معمتعتتمه همه وطمعصهوو رمصهطصم‌وم ممتاهتتامع . 00مصقطصه 6 1۳۱۵ 0عوهوا-مم1/ه 
2022(۰) 

-طعنط ۵۶ وتورامصه معصممصمموم همه ممتامع۳20 * ر(2022) رانا ۷۰ عک طها ین ولا۷۷ .۷ و1۵9 26 بنا۷ ۷۲ بط۷1 .18[0] 
0 واوحتمهماه]۱۷ اعد رومامصون فبامزیه۷ اه ومتهنامج مالدممصومی لهامصصمصصرآمم متمما لهعتتصنانه عمصقم]وه 
۰ ,(11) 31 ,وعتاا5000 

حمصونا0ظ عتصما تم و0مطاه۱ ممتاحتاع مصوتومموامو مه مصلوجمو و ه1لاصزمعو۲۱۵ ,۲۱ .9 4صه صبت نه م۱۷۵۴60152 19[۷۷۰] 
۰ ,3967 ,18 .20 ,19 ,۷۵۱ ,5۵7۱5076 ۲۹۵۷1۵۷۱۷۲۰ ۸ :1۳۱/۵(۰) دومحم ۱۷۲۵2۵1 

0۵۷0۲ اا-عتمنعقا مصنود عمامنتمع ۱۳۸/۵ ۵۶ صمتاهم)ط10۵ تمد مطه‌ظ. پابفطوها .۱ 4صه م2جصصع2 ,20[۳1۲] 
0 1۳0۹۵۵09 ای .عمصهص‌هاتهم تتقطا من قاععگه بواتتصصط مج متنااجمم‌جوها ‏ تمصع ۵ طمتاجهلتقصمم 
۰ ,2723-2732 و11 ۱0۰ ,67 ۷۵۱ ۱۷6۵6۱۵۱۱6 0۳00 ۱5۳۱۵۱۵۵0 

عصتامتمنی »م1 عمج فص لز۵امظ لحتمه کممتننهع ۸ ؟ رزامعمم]۱ ی ک رحع0۵0 بع رتاهطع۱2 .۲ ,2207240 .21[۳1۲] 
۵ دنرگ م۵ ۵6 ۵ ۵۱۵ "روتماهنتامه ماتوومجصمی اهامصصتمصصراو منصمز گم تمتقطهه ‏ عصت‌صهه 


51۲۱/6۲۷۲, ۷۵۱۰ 31, 80. 17, 00۰ 1973-1985, 00۰ 


سال سی و چهار» شماره سهء ۱۴۰۲ نشریه مهندسی متالورژی و مواد 


۹۲ مدل‌سازی سه‌بعدی اثر خستگی در کامپوزیت‌های پلیمری فلزی یونی 


-حمط عم مهمتمرمره 0۷۵۲ امه لگنج مجح متعم1 وممنظ 0مصامامجمع ۸ تمصع .۲ را2۷طعع/.۱۲ ,2۵0720 ,22[۳1] 
۰ ,612283 ,10.4 ,35 ,۷۵1 یرگ ۳26 رطم016162مص ولوعتماووط طذ وتمتمیتامج ۵۶۱۳۸۲۵ ومتاعع1)تاصم10 تهممز 

اعد 501 ۶ ولوولممد اممصه۱۷۲۵۷ عصتصوظ ؟ بصمتمامطهو. .وک 24صصهم . رتصتطه:۲۱۳۵ هر بتعتتها م۸2 .23[۳۴] 
۰ ۵ ۳۲۵۵۲۱۵ م0صانهن - رامصقط) ۴ -عصممتان ۳16 ۲۵۵۲162/60 

مصجماصمصه وتمعه االهمتصمطهمتمماه اع۲مص هچ ۵۶ پمممممام۱۵۲ * رلعمع۱ .و کی مطمصمله بظ رتتاممتطعع820 بک24[1] 
۰ 00۰1/2۰ ,2 .20 و1 ۷۵۱۰ ,5۳۵۵۵۲۵5 ۵۳0 قآ۵ ۷۵۱۵ 5/۱۵۳۲ "رصح /801ع5 طمت1012 ۷ 92800 92)01121 صرق 2۳0 

,1 ۳۵۲۸ :ذابتای متاععای ها ومالوممصصمی مصتصتتهای مصع‌مامصممه مه روامتمعهتامن ۵۶ فصت‌صعظ؟ هتسباو بکد عک. عقوم > [25] 
0۰ ,00.186-198 ,1-2 ,۷۵۱ ,480 ,۷۵۱ ,تاکتصجهطن [جم1) وامصهمتامعاظ ۵۶ امصامز . عمتاممل۱ ۲۵50085 

,80.10 ,104 ,۷۵۱ ,ییا ۸۵/۵0 ]۵ هل روماتدمم‌جومی قاعصه تمصطنرآومر منصما صا فمتصصهصرول معتقطان ,۳۵:۲۲ 26[۳] 
0۰ ,15 1049 

عاتوموصصم آهاع]۱--مصوامظ منصما مه امه 0عوظ-عمنووطاط مجح 4مممزن0امتاصمن) ۸ مصع] ...کت مطقطان ..27[2] 
۰ ,19-529 00.5 و5 .20 ,13 ۷۵۱۰ ۷۵6۱۵0۳6۵ ۵۱ ق0عم۵ی 1۳۵7 ۲۲۱۲/۵۵ ۰ 200۲9اهظ 

تعصصرامم متصمد ما قممتامماه همه مموتاو حاونا۲. روما .۲ .1 عک رقاعلظ .۲ .ظ رتعع۵۲م5تمصصااه ۷۷ .۲ رنط۲۱۵06 ۷۷۰ ,۸ .1 .28[۳] 
۸۵ ۶ اوصتاه؟ رماتعطهها لهعتصصمطم-متاممام فطا که فصمتاهعتامع1۳۲ لماجمنموون مه لهعزتمصصتاه. عتممنل‌قمهه 
1۰ ,074509 00۰ ,80۰7 ,109 ۷۵1۰ ,۲۲۷۹169 

۲ 200۳020 ۱6۷۵۲۱ لحتامط اهلم)تانج 40ج منعم1 و22 0مصاصاصصمی خر ملظ ,۲ مه رفطمعلا .۲ ,22724 .29[۳] 
,.00.612283 4 ,0ظ ,35 ,۷۵۱ کرک 6۲ 01۴162010 محر فتقعتماووط ط عتماعتتمع ۲۱۳۸/۵ ۵۶ طمتا1060416162 عومصتلممه 
2018۰ 

) مهد جع رلهتاممامم لقع :ومام‌عهد اقامتانه قج ومتقنطامه (رطمظ) تمصصراوو وتامهمتا0عاظ۳۳ رصمطمن-ت2ظ 30[۷۰] 
٩۳۱۳ 0۲659,‏ :۷۷۸۵ متطقطعصتااهظ .(1-765 .00 ,136 

۲ اوعد عالوومرجصیی اقامصه راوج متصما رم حمتاعقصمم‌طصمی 010017۵عت 00وجها 061مص ۸۵ مضع] ۷۰ 4ص2 رععه(1 ب31[۳۴] 
۰ ,00.43-49 ,224 ۷۵۱۰ رامع« ۸۰ معط ۵۳۱0 عمعوعی * رتصمجصهتاعحعطط کصمجمرعمه مرول 

7 .۰ 2001102110105 ۵01621 4ج هرفن1۳ ۲۱۷۰ :ومازومم‌جمی آهامصمصوراوم منم متا .۲ مک 24 تمممرصنطقطه .32[3/1] 
,۰ ,00.197 ,1 ,0ظ ب14 ۷۵۱۰ ,5۱۳۱۵۵۲۵۲۵5 ۵۳۱۵0 دا۵ز111016۳ 

۲ 5099018 مالومم‌ممه هام تمصصرامی منعمد عم ممتامقممم‌صصرمی مبتامن0۳۵0 2800 عم ۴۸۵ ص12 ۷۰ کی قوهمنا بک[33] 
۰ ,43-419 00۰ ,224 ۷۵۱۰ ,60و ۰ ماگ له کرمعوهی ۰ اصمجه تباعععمطط کصمجصطممه۱ مرول 

تصمتامعصوا اهعتتامعا0مصمطمعصه مالوممصمی لمامصص‌مصصرامم متصما * و جنک .لیکا روماماقظ هو عصتاگ ‏ یک و ۳۵۵2 ,«34[1] 
0۰ .79-289 00.32 ,20.3 59 ,۷۵۱ ۲۱ ۵۲ روع615060۷ج 4طع ۲۵۷1۵۷۷ 

مج کعماهتامعه مالومموصم امامص-مصصرآهوم متما ط وعمصع207 کصعمم ؟ . هعلهفه یک ر تک ,لب یکا-1 نام بقع .1 ,10 35[)0] 
۰ .00۰1037-1066 ,7 .20 ,38 ,۷۵۱ ,66 ۵ و۳۳۵ رعطمتاه 2001 40ص2 عطتاملمصه تتفط 

2 66 ۸۱۱01 .06 , مصلهناگ .1 , 216و بک , ۱۷۱2212 هم صتاوع۲ .۲ , 50وماظ ب), [۷۱۵2 .۳۲ , وعصطلجظ ۷۰ , نک کار عصتصصنظ نه(36] 
4:6913)2014(۰ وم 901 روصمتاهنآممه ممهمر6 ص ومامفباصط أهلم)تایه تمصراوم و۷ تام2متامعاه متمم[ > 

عصتلصعه موعقا ما ممتای‌هاه: م02 ۱۳۸۲۵ مصتمتهتاومط. ,۷ مصمهاه .از مه رلدزع-01صموه]۱۷ ,۷۲ بتصهاوع‌طمصصم ۳71[ 37] 
1۳0۱5۵61015 ۳ .01600۳0065 0مصماادن رها ممصحماساونل صهتوعنیهن 10621 1م6002ع]-ممظ عمتراممه :عاجمجطممهامفت 


۳2661۲۵0۷ ۵۷6۵6, ۷۵۱, 63, 20. 4, 00۰1689-1695, ۰ 


[۳۸ م محمودنژاد. ۸ هفته چهارم دوره جامع آموزش عمومی کامسول- مش بندی در کامسول. 


نشریه مهندسی متالورژی و مواد سال سی و چهار» شماره سه. ۱۴۰۲ 


آرزو خراجی, ادیا نقوی, حجت زامیاد ۹۳ 
۵۷۵ ]۱۷-ع۲منوها مصتوا وعمتعتامج ۱۳۸ ۵۶ ممتاهم)م10۵ تمرمدهظ . مدا ۲ ی مم2جصصه2 .39[۳] 


۱ 1۳050605 لتایل "رومصهص‌همتتهم تتقطا هم قاععره بواتلتصنط مج متبااجمم‌جها ‏ تمصع ۵ ممتامهل8صهم 


۰ ,2723-2732 00۰ م11 0۰ ,67 ,۷۵۱ ۷۵۵۵۱۵ ۵0 ۱۱6۱۳۱۵۱۱۵۵0 


سال سی و چهار» شماره سه, ۱۳۰۲ نشریه مهندسی متالورژی و مواد 


۹۴ مدل‌سازی سه‌بعدی اثر خحستگی در کامپوزیت‌های پلیمری فلزی پونی 


نشریه مهندسی متالورژی و مواد سال سی و چهار» شماره سه, ۱۴۰۲ 


